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Введение 


Актуальность задачи. При разработке месторождений полезных ископаемых, в 
частности, угольных пластов, важное значение имеет задача снижения интенсивности 
проявления негативных свойств разрабатываемого массива. Одним из основных спо- 
собов решения задачи является целенаправленное изменение состояния массива путем 
внешнего воздействия. В условиях анизотропии обрабатываемой сплошной среды ка- 
чественная реализация процесса обеспечивается введением автоматического управле- 
ния параметрами воздействия, при этом состояние массива представляется с помощью 
математической модели. Неопределенность характеристик среды приводит к необходи- 
мости использования элементов искусственного интеллекта (в частности, элементов не- 
четкой логики). В этой связи тема работы является актуальной. 

Цель работы — обоснование структуры математической модели и алгоритма ее 
реализации. 


Математическая модель нелинейно-упругого режима 
фильтрации жидкости в трещиновато-пористой 
сплошной среде 


Г идравлическая обработка угольного пласта осуществляется путем напорного на- 
гнетания жидкости через скважины, пробуренные по пласту согласно технологическим 
схемам, принятым в соответствии с горно-геологическими условиями. 

Исходная модель строится на основании уравнений движения и неразрывности 
ДЛЯ ЖИДКОСТИ ПУ, в которые входят скорости движения жидкости и смещений твердых 
частиц: 


ОР 

е- тб -и,); (1) 
бт 9 
+ Им =0. 2 
а ® © 


В сцементированных средах, к которым относятся угольные пласты и вмещающие 
породы, можно пренебречь скоростью перемещения твердых частиц по сравнению со 
скоростью жидкости всюду, кроме уравнений неразрывности. На основании резуль- 
татов исследования этих уравнений в режиме фильтрации жидкости с относительно не- 
большим (до 103 ат) перепадом давлений в работах [1], [2] получено уравнение не- 
линейно-упругого режима фильтрации: 


ОР д ОР 
—=и— 1 +а(Р-Р)|—Ъ 3 
я 5 [-+0(Р-Р,)| ах (3) 
А 
где у = о коэффициент пьезопроводности; 


№ "а, +а„) 


а=а, +а, ар, 


а,›а,,а„›а, — соответственно коэффициенты изменения проницаемости, пори- 


и? 
СТОСТИ, ВЯЗКОСТИ И ПЛОТНОСТИ. 
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Полная математическая модель процесса должна содержать комплекс матема- 
тических зависимостей, на основе которых можно представить ход процесса и рассчи- 
тать его параметры. В данном случае для получения математической модели процесса 
гидродинамического воздействия на пласт уравнение (3) должно быть дополнено на- 
чальными и граничными условиями, которые, в свою очередь, определяются техно- 
логическими схемами и режимами обработки [3], [4]. Значит, постановка задачи 
исследования и расчета параметров процесса формируется в соответствии с конк- 
ретной технологической схемой. 


Основные варианты постановки задачи моделирования 


Предварительно необходимо сделать несколько замечаний по поводу терми- 
нологии. При воздействии на пласт через группу скважин говорят о галерее, батарее 
или каскаде скважин. Некоторые авторы применяют тот или иной термин в зависи- 
мости от схемы расположения скважин, другие же оперируют этими понятиями 
произвольно. 

Предлагается при обработке угольных массивов связать эту терминологию с тех- 
нологией воздействия и далее говорить о галерее скважин, если нагнетание произво- 
дится поочередно из различных скважин данной группы независимо от их взаимного 
расположения, под каскадом же понимать одновременное воздействие на пласт через 
группу скважин в условиях гидродинамического взаимодействия потоков, обеспечи- 
вающего равномерное насыщение массива. При такой терминологии задача расчета га- 
лереи не отличается от задачи для одиночной скважины, при каскаде же имеются 
принципиальные отличия, связанные с интерференцией потоков и взаимодействием 
контуров. 

При постановке задачи исследования для различных схем уравнение (3) допол- 
няется соответствующими начальными и граничными условиями. 

В начальный момент времени / = 0 задается распределение давления в области 
фильтрации: 

Рен] =Р.. 4 
На скважине задается режим давления 
Вы 


или темп нагнетания 


Е ОР) = са (и) (5) 
т 


где Рс(@) — давление на скважине; х; — координаты скважины; Сс — контур 
скважины; 4(1) — темп нагнетания; с — коэффициент, зависящий от размерности. 

На границе области фильтрации задается либо давление, если область оконту- 

рена скважинами или выработками, либо условия непроницаемости на контакте с 

боковыми породами или на «бесконечности» (когда обрабатывается нетронутый 


|@р 


массив): 
Р(х, 1], =Р, (6) 
или 
ОР 
0, @ 
ди 


где /'— граница области фильтрации; Рг- давление на границе области. 
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Таким образом, сформулированные задачи исследования параметров процесса 
гидродинамического воздействия на пласт исчерпывают основные рациональные схемы 
обработки. Технологические схемы с иным расположением скважин будут являться 
частными случаями рассматриваемых задач. Дальнейшая работа заключается в матема- 
тической формулировке каждой из отмеченных задач и решении полученных уравнений. 


Моделирование процесса воздействия на пласт 
через одиночную скважину, пробуренную впереди забоя 


На рис. 1 приведено схематическое изображение варианта воздействия на пласт 
через одиночную скважину, параллельную линии забоя, в режиме фильтрации. 

В этом случае скважину можно считать расположенной в начале координат, 
другие варианты положения скважины можно свести к этому. 

Основу модели составляет уравнение (3). 

Начальные и граничные условия формируются из следующих соображений. При 
данном расположении скважины обрабатываемая зона ограничена, т.к. при длительном 
нагнетании жидкость выйдет на выработки, оконтуривающие массив угля, подвергаемый 
воздействию через скважину. Поэтому область фильтрации по оси ОХ ограничена от- 
резком ОГ., причем величина Г, может либо задаваться (когда известны размеры зоны, 
которую необходимо обработать одной скважиной), либо рассчитываться (когда задано 
ограниченное время обработки Тобр). 


Рисунок 1 - К расчету параметров воздействия на пласт через скважины, 
пробуренные впереди очистного забоя: 
а) одиночная скважина; 6) каскад 


Следовательно, в качестве начальных условий задаются: 
а) давление на скважине в начальный момент 


РИ], = Р(0,0) = Р.; (8) 


6) распределение давления на отрезке (0,Т] при #=0: 
Рок) ‚„ =Р(>,0}; (9) 


По смыслу, начальное распределение давления, очевидно, равно давлению 
газа в пласте: 
Р(х,0) =ВБ; 
Краевые условия. 
На левом конце (х=0) задается давление на скважине или темп нагнетания, 
причем эти величины могут быть либо постоянными, либо являться функциями вре- 
мени: 


Р(х,г) 


хо = Р(0,1) = 2. (1); (10) 
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ИЛИ 


а Е (р 
х=0 

где 4 - удельный расход (расход на единицу площади фильтрующей поверх- 
ности); с — коэффициент приведения, зависящий от размерности величин, входящих 
в выражение (11). 

На правом конце (х = Г) условия определяются схемой, а именно: 

а) если на расстоянии Г от нагнетательной скважины пробурена отточная и 
исследуется движение жидкости от нагнетательной к отточной скважине, то в точке 
х=[ задается постоянное давление, равное давлению на выработке (обычно считают 
его равным атмосферному): 


Р(х,г) 


р = РИ =Р,. (12) 
[> 


Такое же условие задается, если рассматривать фильтрацию в сторону выра- 
ботки, но при этом надо учитывать влияние опорного давления [2], [3]; 

6) если исследовать движение жидкости по мощности пласта в сторону кров- 
ли (почвы), то на границе задается условие непроницаемости: 


—- = 0. (13) 
И 

Решая уравнение (3) с начальными условиями, выбранными из (8) — (13) в со- 
ответствии с конкретной задачей исследования, можно найти распределение давления 
на отрезке [0,Г.] в любой момент времени, что дает возможность рассчитывать пара- 
метры соответствующей схемы воздействия на пласт. При этом, как указывается в [2], 
решив задачу при некоторых нормированных условиях, можно получить решение боль- 
шого класса задач, переход к которым осуществляется обратным пересчетом нормиро- 
ванных решений. 

Уравнение (3) — это нелинейное параболическое уравнение, точное решение ко- 
торого весьма затруднено. К настоящему времени разработан ряд численных и анали- 
тических методов для приближенного решения уравнений этого типа. Большой вклад в 
разработку и практическое использование уравнений этих методов внесли ученые 
Института кибернетики с ВЦ АН Узбекистана под руководством проф. Ф.Б. Абута- 
лиева. Разработанные ими методы и программы успешно используются при исследова- 
нии процессов в области мелиорации, добычи нефти и газа и т.д. [4]. Эти методы в даль- 
нейшем используются в работе. Однако различия в условиях воздействия, целях иссле- 
дований (в частности, разработка инженерной методики расчета параметров обработки 
угольных пластов) определяют необходимость проведения самостоятельных исследо- 
ваний для угольных пластов. 

Как отмечается в [2], проницаемость пласта различается не только вдоль различ- 
ных осей. Структура пласта весьма сложна, и картины фильтрации, полученные для 
постоянных усредненных значений фильтрационных параметров, дают идеализирован- 
ное представление о характере процесса. При такой форме области воздействия распро- 
странение влаги должно происходить равномерно, чего на самом деле не наблюдается. 
Следовательно, необходимо так подобрать параметры модели, чтобы получить реаль- 
ную картину процесса. 

Один из путей подбора таких параметров — задание анизотропии пласта по 
всей области фильтрации независимо от направления осей, при этом, естественно, 
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необходимо учесть, что коэффициент проницаемости вдоль оси 07 в любом случае в 
несколько раз меньше проницаемости по простиранию (оси 0х). Поскольку пласт 
состоит из блоков с разной проницаемостью, пересечен трещинами различного про- 
исхождения, то в нем можно выделить некоторые области с различной величиной 
коэффициента фильтрации, и соответствующим образом задать набор коэффициентов 
проницаемости. Ввиду того, что выделить на практике такие области, не проникнув в 
пласт, невозможно, в первом приближении можно задать набор коэффициентов 
проницаемости хаотично, используя таблицу случайных чисел. Благодаря этому приему 
можно будет исследовать процесс фильтрации в условиях, более близких к реальным по 
сравнению со случаем, когда проницаемость осредняется по трем (двум) направлениям. 

На рис. 2 приведены результаты решения задачи в виде рассчитанных линий 
равного давления (в процентах от давления на скважине) в области воздействия при 
постоянном значении (рис. 2а) и случайном наборе значений коэффициента проница- 
емости (рис. 26), т.е. коэффициент проницаемости представляется переменной, имеющей 
стохастический характер. Как видно из результатов моделирования, жидкость распро- 
страняется по всему пласту неравномерно, образуются слабо увлажненные участки. При 
этом жидкость может найти путь фильтрации, по которому быстро выйдет на соседнюю 
скважину, оставив большие блоки неувлажненными. Обходя эти блоки, жидкость «за- 
купоривает» имеющийся в них метан, а это может привести к образованию зон с вы- 
соким содержанием газа под давлением, и увлажнение даст нежелательный эффект в 
отношении внезапного выброса. 

Таким образом, используя технологию нагнетания жидкостей через одиноч- 
ные скважины, добиться равномерной, качественной обработки массива практически 
невозможно. 


—= — > 
а) Одиночная скважина 6) Одиночная скважина, 
К=Сопя К=Гаг 


Рисунок 2 — Результаты моделирования процесса воздействия через одиночную 
скважину при постоянном и варьируемом значениях проницаемости 


Выводы 


Процесс гидровоздействия на сплошной трещиновато-пористый угольно-по- 
родный массив адекватно описывается нелинейным уравнением параболического типа, 
решение которого производится с помощью метода конечных разностей. 

Для моделирования процесса воздействия в условиях неопределенности зна- 
чений параметров объекта обработки предложено ввести характеристику проницае- 
мости среды в виде стохастической переменной, что позволяет получить реальную 
картину процесса и эффективно использовать результаты моделирования при по- 
строении алгоритма управления технологией гидровоздействия. 
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И.М. Рашузй, А1-Лег Е Огъа 
Тйе АррПсаноп о} Матетайса! Мо4аейие т !Ше Зуяет оГАиютайс 
Сотто/[ о} Ргосез5 ор Нуагаийс Асноп оп СтасК-Рогоих Мазяме 


т Сопашопз о таейийе Рагатаету ОГбиггоип те$ 

ТБе ргоМет оЁ сопзбласНоп оЁ зузет оЁ сопёо[ оЁ 1есппоюозлса| ргосез$ оЁ Вудгаийс 
асбоп оп ипеуеп газ сощашише соа| зеат аз а \У’ау оЁ диесеа уапайоп оЁ Из зава$ 1$ 
сопз14егеа. ТВе таетайса| по4е|, Базе оп попПпеаг едиайоп оЁ рагабо1а фуре, ргорозе4 ог 
дезспшо оЁ ргосезз. ТВе еззепна! Ееавге оЁ ргосезз 15 шдеНпце оЁ рагатеегз, ай гергезепЕ 
Бу гапдот уапаЫез. 

Сотр!ех тафетайса| то4е| плиз сопашз Фе зузет оЁ тафетайс$ дерепдепб, 
УЛись аге Фе Базе оЁ гергезетаноп оЁ 5айа$ оЁ ргосез$ ап сайсиае Из рагатеегз. ш ргезеп 
сазе Юг оМашше оР тафетайса| по4е| оЁ ргосезз оЁ Вудго-4упаписа| асйоп оп зеат Фе 
Баз1с ефиайоп ти Бе а44еа Бу ша! апа Боипдагу сопоп$, фа{ 4ейпе ассог@те. ю {ес- 
поостса| зспетез ап фурез оЁ ргосеззшо. Зо, ше Ююипдайоп оЁ Ше ргоет оЁ шуезйеайоп 
ап4 са1сшайоп оЁрагатеегз оЁргосез$ плиз Бе 4еаЙ ассог4то {0 отуеп 1есрпо]остса!| зспете. 

Те соа] зеат Ваз сотрИсае4 збгасвге, сощатз ЫосК$ УЛ уапои$ Но\’ соеЁйает, 
ап4 уюте паи ещегз птеоцаг. Еог рго\1А те оРеНесйуе асНоп \уе Вауе 1ю Нпа ФЧза1БиНоп 
оЁ сБагасетз, Биё И 1$ парозЫе улфоиЕ дезбасйоп оЁРтазяуе. 

п Фе асе Ф1561Бибоп ог Ноулио рагатеег$ оЁ зеат аз а гапфот уапа ез, шпйеа Бу 
тахипа! ап тпшила! уа[чез, 1$ ргорозед. ТЫ$ тео ргоу14ез $ ю шуезНоае Ше ргосез$ аз 
теа| оБ]ес{ ап4 гу ю Вп4 а угауз Юг оуегсоте оРпесайуе шРаз1оп оЁ 1$ ргорейу оЁ соа! ап4 
тоск таззуе. 

Те гези№$ оЁ шуезйеаНоп$ зВоу’, Фа изше фесбпоо1са| зсВетез оР шЙмепсе у 
шее сСышК$ итоги диаШануе ргосеззто 15 ппрозз1Ые. 

Тье угау оРа п5е оЁ еНесйуепезз оЁ ргосеззшо 15 шиодасНоп оЁ4есбпо|огу УЛ этоир 
оРсЬК$ рго\14е4 соп4оп$ оЁ ище[егепсе оРсопёгагу збеапл$. 


Статья поступила в редакцию 08.12.2012. 
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